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The present invention relates to a method for 
determining position and velocity of targets from 
signals scattered by the targets. One uses a first 
and a second station comprising 
transmitter/receiver (si ,s2) of electromagnetic or 
acoustic signals and carries out mono-static 
measurements (M1,M2) from each station, and 
also a bi-static measurement between the 
stations (812), and calculates first through the 
two mono-static measurements a number of 
target candidates with 2-dimensional position and 
2-dimensional velocity, whereupon these are 
tested against the result of the bi-static 
measurement and the target candidates which 
are found In all measurements with suitable error 
margins are retained. 
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Method for Determining Position and Velodty of Targets from Signals scattered by 
the Targets ; • 

The present Invention relates to a method for using signals scattered by one or more 
5 targets for detemilning position and velocity of each target. The invention is 

especially useful In the case where the positions in space and velocity vectors for a 
very large number of targets are to be determined. The invention is' for instance 
intended for use with the system for detemining position and velocity of targets for 
which this day a patent has been applied with Hans Hellsten as inventor, and which 
10 patent application 0101661-7 Is hereby incorporated by reference when regarding a 
more detailed reasoning of the system of transmitter and receiver, henceforth called 
stations. 

The invention is based on a 3-dimenslonal position space and a 3-dimensional 

15 velocity space, which together form a 6-dlmensional state space with up to ~-10^ 

cells, where in principle each cell shall be assigned the value 0 (no target) or 1 
(target). The large number of cells In the state space place special requirements on 
the signal-processing method, and one possible method is dealt with In said patent 
application by Hellsten. In the present patent application an alternative method is 
20 suggested. 

Due to the large number, ~ 10^° . of cells in the state space, known projection 
methods are slow and cumbersome. The simplest known method is where all 
conceivable target positions are formed for three stations and N targets. They will 
25 be fewer than in number, and it is easy to detemilne position as well as velocity 
vector for the target for every possible target position. In this way - candidates 
are obtained, each of which is then supported or discarded by further independent 
measurements. This gives a method of processing which requires K N' operations. 
If K Is large then this method can be much too slow for Interesting values of N. 

30 

A specific property of the systems for which the use of the present method is 
Intended Is. according to the above, that the information In an Introductory, detecting 
phase can be regarded as binary, which can be utilised. Further the Input Is from 
the beginning distributed between different stations. This means that a Ideal, 
35 distributed processing which calculates parameters of position and velocity for the 
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targets gives an advantage, as one then receives a paralleling effect on the sensor- 
near level. A comprehensive communication between stations is then necessary. 

The stumbling block is the association problem, which is to in a correct way pair 
5 together data from different radar stations. An association where it is unclear 
whether it is correct or not can be called a candidate. Gradually it becomes dear 
whether a candidate is a correct association, that is a target, or a false association, 
sometimes called a ghost. The problem is then to handle the initially very large 
number of ghosts, that is to find the targets among all candidates. 

10 

The purpose of the invention is to solve this problem, which is done by the invention 
being given the features that appear from the following independent claim. Suitable 
embodiments of the invention will be evident from the remaining claims. 

15 The invention will be described- in more detail -in the following with reference to 
enclosed drawings, where 

fig. 1 shows how the geometry around two stations on a plane and how 
rotation-symmetrical conceivable target positions on measurement 
20 from two stations can be transfen-ed to points on a half-plane, 

fig. 2 shows how velocities in a 3-dimensional space can be transfen-ed to 
non-unique velocities in said half-plane and 

fig. 3 shows how one can determine position and velocity by combining 

measurements from two pairs of stations that have a common station. 

25 

If one studies figure 1 with two stations, St, Sa, on the ground, which see one and the 
same target, f,, in mono-static measurements, one realizes immediately that they 
place targets along a half-circle in the position space above the ground. Bl-static 
measurements between the two stations share this symmetry, so that a half-circle 

30 can be accepted or discarded in its entirety. For each point on the half-circle, an 
estimated velocity v is detennined up to a line In the velocity space. That Is, two 
stations together ascribe a candidate (target or ghost) a two-surface In the state 
space. A famfly of such two-surfaces in the state space fonn a relatively compli- 
cated sub-set In the state space. It is therefore desirable with an Introductory 

35 calculation where candidates can be represented by points in a reduced state space, 
rather than by areas. In the case where one studies position as well as velocity, as 



wo 02/093191 



PCT/SE02/00894 



10 



15 



described above, the reduced state space will consequently be 4-dimensiona!. 
Further some type of candidate reduction should occur in this reduced state space, 
so that the final processing In the full state space becomes simplified. 

The basic idea of the invention is to work in low-dimensional. < 6. state spaces that 
are common to several sensors. Optimally the areas in the full state space, corre- 
sponding to certain sensor registrations, will be represented by points In the lower- 
dimensional ditto. As the reduced state space is common to several sensors, one 
can under certain conditions obtain a sensor-near and local elimination of many 
ghosts by calculations which are not demanding to process. Such a reduced set of 
candidates, together with other similarly reduced sets of candidates make up input 
to a final and complete association in the full state space. This final calculation is 
then, by the relatively small size of Input amounts, also sunnountable in ternis of 
processing Intensity. 



From the symmetry in figures 1 and 2 one can see that one can allow the points on 
the arc corresponding to a candidate f, be represented by a point P(ti) In a half- 
plane P,2 on one side of the symmetry axis. In figure 2 one can see how the velocity 
V can be determined from mono-static and/or bl-static measurements from a pair of 

20 stations by Doppler measurements, except for a component in tangent to the arc 
direction. This means that the velocity lies on the line Lrs,.s2.f,.v;. All velocities on 
this line will be represented by the same 2-dimensional velocity u. which is a vector 
In the half-plane P„. In reverse, it applies that if the velocity In the half-plane P„ is 
known, one knows that the velocity lies on the line L(si,S2.ti.u). where u is now 

25 regarded as the canonic representative of u in three dimensions. 

In contrast to the simple problem of finding the common line between measurements 
from two stations, tiiat is the arc from the intersection of two (half-)spheres. or the 
common point between measurements from three stations, that Is the point of inter- 
30 section of three (half-)spheres, the common points from two or three bl-static 
measurements between the same stations are more complicated to find, as each 
measurement generates a (half-Jellipsold with the stations at the focal points, it Is 
more difficult to calculate common points between two ellipsoids Uian between 
spheres. 

35 



wo 02/093191 ^ PCT/SE02/00894 

4 

The Invention utilises therefore mono-static measurements between pairs of stations 
to calculate candidates. It then uses the associated bl-static measurement to 
eliminate the ghosts. In this way one does not need to calculate common points 
based on the ellipsoids, but rather loolcs only at the points where the mono-static 
5 measurements indicate that there are candidates and sees whether the bl-static 
measurements have one or several candidates at the distance in question. This is 
much simpler than calculating common points fi-om bl-static measurements. In 
cases where the distance-based analysis does not give sufficient elimination of the 
number of ghosts, one can use as well the compiled Doppler information. The 
10 Doppler measurements share acconling to the above the symmetry of the distance 
measurements and by requiring compatible mono-static and bl-static Doppler 
measurements a further reduction in the number of ghosts is obtained. 

In the following a more detailed description is given of the method in connection with 
15 an amount of distributed Isotropic transmitters/receivers of electrorha^^^ 

acoustic signals. In the concrete example we consider a number of radar stations. 
Each radar station sends out a signal that is received by the station itself and>also by 
sunx)unding neighboura. We refer again to the above-mentioned patent application 
by Hellsten for a more detailed description of how the system can be embodied in 
20 general. 

In a first step one uses a first s, and a second 52 transmitter/receiver of electromag- 
netic or acoustic signals and carries out mono-statlc measurements iW, and Mafrom 
each transmitter/receiver, and also a bl-static measurement B,2 between the 
25 stations. First the two mono-static measurements create a number of candidates 
with 2-dimenslonal position and 2-dimenslonaI velocity. Then these are tested 
against the result of the bl-static measurement and the candidates found in all of the 
measurements with suitable emor margins are retained. 

30 In a second step one carries out corresponding calculations with one of the two 
original stations, say s,. and a third station s,. The same calculations that were 
carried out for s, and 82 are carried out for s, and S3. One then sorts, for each 
calculation, conceivable targets according to the distance to the common 
transmitter/receiver, here Si, and determines the distance bins In both calculations 

35 that simultaneously contain one or several candidates. For these candidates the 
radial Doppler velocity in the two measurements Is studied. If the difference is 
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smaller than a predetem^fned value then the candidates are retained; if not then the 
are discarded. • 

Rgure 3 Shows how a target is represented by a point In the half-plane P„ and 
another in the half plane P,,. The 3<Iimensional position is given by the point of 
.ntersection of the half-circles. Thetrue velocity ^ for Ms represented by the 
velocity uu which is a vector in the plane tangent to the point P(tO in the half-plane 
P„ and corresponding vector u, in the half-plane P^r. Each such velocity can be 
represented according to the above on a line In the 3Klimensional space The true 
velocity V lies in the intersection between these lines L(s,s,t,u,) and L(s,s,tuU^. 

In a third step one can carry out a third bl-static measurement B,, between the two 
transmitters/receivers, here and s, which is not the common transmitter/receiver 
according to the above. The result can be used together with the candidates that 
remain after step two above to- further reduce. Wnumber of ghosts 

The local treatment around three stations has through this procedure, with initially 
paired calculations around two stations and then Including an three stations in the 
calculations. In an elegant way requiring relatively few calculations, given a number 
of targets and possibly, as a consequence of imperfections In the measurement 
system and the application of it. some remaining ghosts. It is interesting that in 
parallel with this, in a larger system, a number of other stations can also have given 
a number Of targets. The parallel handling is very favourable and an absolute pre- 
requisite according to the above for the processing to be at all possible in a larger 
25 system with many Stations. 

To improve the accuracy In a,e calculation, and discard pcssWe remaining ghosts 
one can of course allow further stations. v*lch have a sufficient range to be 
induded m the calculations. On. can consider further n«,no-statlo a, well as further 
Di-static calculations. 

ft ta. almady been mentioned several flme, that the p,esent InvenUon Is Intended to 
be Included In a brger .jrstem of station,, for example that which Is presented In the 
saw patent appOcafion b, Hans HeDsten. In this application H b stated that a 
number of station, can be placed as gdd point. In an essemially equidistant grtd on 
a .urface. ««ch «mlU the su,ve,«l «ate space, e.,. a g™und surface. Tl» Wea In 



20 
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the patent application of Hellsten is that the ranges of the stations shall be such that 
each target is detected on bi-static measurement by at least six independent sets of 
transmitters and receivers. For this reason the range of the stations is 2cf, if the 
distance between the stations is d, if we consider an essentially planar surface. If 
5 the surface is not essentially planar then the ranges must be adapted so that an at 
least 6-fold bi-static overlap of frie surveyed state space is nevertheless achieved. 
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Claims: 

1. Method for determining position and vefocity of targets from signals scattered by 
the targets, c fi a r a c t e r i s e cf In that one uses a first and a second station 

5 comprising transmitter/receiver (si.Sz) of electromagnetic or acoustic signals and 
carries out mono-static measurements (Mi.Ma) from each station, and also a bi-static 
measurement between the stations (B12), and first from the two mono-static 
measurements calculates a number of target candidates with 2-dimensional position 
and 2-dimensional velocity, whereupon these are tested against the result of the bi- 

10 static measurement and the target candidates found in all the measurements with 
suitable en-or margins are retained, whereby the calculations take place in a 2- 
dimensional position space and a 2-diinenslonal velocity space and give the position 
and the velocity modulo rotational symmetry around a line through the two stations 
(S1.S2) In that one utilises that rotational symmetry results in the ability to represent 

15 - conceivable measured target positions In a 3-dimensional position space, which lies 
along arcs with centre on the axis of symmetry, by a point {P(f)) in a half-plane (P«) 
on one side of the axis of symmetry, which makes a 2-dimensional position space, 
and that the velocity can be represented by a velocity In a 2-dimensional velocity 
space that Is a plane tangent to said point. 

20 

2. Method according to claim 1 . c /i a r a c t e r / s e c/ in that one establishes the 
position for a target candidate in a complete 3-dimensional position space and 3- 
dimenslonal velocity space through con^sponding calculations with either the first or 
the second station (Si.Sz) and a third conresponding station (S3), whereupon for each 

25 calculation one sorts target candidates according to distance to the common station 
and detennines for the distances In both calculations which both contain one or 
several conceivable targets, whether the candidate shall be retained as target or 
discarded depending on the agreement of the radial velocity. 



30 



3. Method according to claim 2. c /i a r a c t e r / s e d In that one cames out a 
third bi-static measurement (Ba) between the two stations (sj.Sj) which are not the 
common station according to claim 2 and retains a target candidate if agreement is 
achieved, othenvlse the target candidate is discarded. 
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4. Method according to claim 2 or 3, characterised \n tliat one uses 
additional stations for mono-static and/or bi-static measurements which are 
compared with the result of the earlier ones. 

5 5. System according to anyone of claims ^-4, characterised in that said 
stations (si.sa.sa) are placed as grid points in an essentially equidistant grid on a 
surface, vWiich limits the sun/eyed position space, e.g. a ground surface, with the 
distances between the grid points essentially the same size, d, and where the range 
of the signals when having an essentially planar surface are at least 2d, implying at 
1 0 least 6 independent bl-static configurations per grid point and that the range, in the 
case that the surface is not essentially planar, is adapted to give the same number 
of bi-static configurations as in the planar case. 
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